Exact solvability of the quantum Rabi models within Bogoliubov operators 
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The quantum Rabi model can be solved exactly by the Bargmann transformation from real 
coordinate to complex variable recently [Phys. Rev. Lett. 107, 100401 (2011)]. By the extended 
coherent states, we recover this solution in an alternative simpler and perhaps more physical way 
without uses of any extra conditions, like Bargmann conditions. In the same framework, the two- 
photon Rabi model are solved exactly by extended squeeze states. Transcendental functions have 
been derived with the similar form as those in one-photon model. Both extended coherent states and 
squeeze states are essentially Fock states in the space of the corresponding Bogoliubov operators. 
The present approach could be easily extended to study the exact solvability or integrability of 
various spin-boson systems with multi-level, even multi-mode. 

PACS numbers: 03.65.Ge, 42.50.Pq 42.50.Lc 



I. INTRODUCTION 

Matter-light interaction is fundamental and ubiquitous 
in modern physics ranging from quantum optics, quan- 
tum information science to condensed matter physics. 
The simplest paradigm is a two-level atom (qubit) cou- 
pled to the electromagnetic mode of a cavity (oscillator) , 
which is called Rabi model In the strong coupling 
regime where the coupling strength g/ui (w is the cavity 
frequency) between the atom and the cavity mode ex- 
ceeds the loss rates, the atom and the cavity can repeat- 
edly exchange excitations before coherence is lost. The 
Rabi oscillations can be observed in this strong coupling 
atom-cavity system, which is usually called as cavity 
quantum electrodynamics (QED) [2]. Typically, the cou- 
pling strength in cavity QED reaches g/uj ~ 10~^. It can 
be described by the Jaynes-Cummings (JC) model [Sj] 
where the rotating-wave approximation (RWA) is made 
and analytically closed- form exact solutions are available. 

Recently, for superconducting qubits, a one- 
dimensional (ID) transmission line resonator Q or 
a LC circuit [5|-|8| can play a role of the cavity, which 
is known today as circuit QED. More recently, LC 
resonator inductively coupled to a superconducting 
qubit (9l-[Tl| has been realized experimentally. The 
qubit-resonator coupling has been strengthened to ten 
percentage. In this ultrastrong coupling regime of 
circuit QED, the evidence for the breakdown of the 
RWA has been provided by the transmission spectra . 
The remarkable Bloch-Siegert shift associated with the 
counter-rotating terms also demonstrates the failure of 
the RWA [l0|. So the quantum Rabi model (QRM) has 
been revisited by many authors. 

In the representation of bosonic creation and annihila- 
tion operators in the Bargmann space ,12.1 of analytical 
functions in a complex variable, Braak [l3| recently de- 
rived a transcendental function in the QRM, which is 



defined through power series in the coupling strength 
with coefficients related recursively. Zeros of transcen- 
dental functions can give exact eigenvalues. By a pro- 
posed criterion for quantum integrability, Braak further 
shows that the QRM is integrable due to the parity sym- 
metry. However, the derivations seems to be outlined in 
a mathematical way. It was also suggested [ij] that an 
intense dialogue between mathematics and physics would 
be needed. In other words, it is useful to shed some phys- 
ical insights into Braak's mainly mathematical solutions. 

In this paper, without the use of any extra conditions, 
like analyticity of the eigenfunction in Bargmann repre- 
sentation, we alternatively re-derive the same transcen- 
dental functions as in Ref. [isj quantum mechanically 
within the extended coherent stats [H, Both zero 

bias and biased QRM can be treated simultaneously. Our 
method is more intuitionistic, and may be more easily un- 
derstandable. The key procedures can be described with 
a simple tutorial, and therefore it is straightforward to 
extend to study the exact solvability or integrability of 
various spin-boso n sy stems. The extension to the two- 
photon QRM [13, ITsI is performed in this paper as the 
first example. 



II. THE QRM WITHIN BOGOLIUBOV 
OPERATORS 

A. Re-derivation of Braak's solution 

The Hamiltonian of a generalized QRM can be describe 
as follows 

H = -]^{e<7^ + /^a^) + a'^a + g{a) + a)<7,, (1) 

where e and A are qubit static bias and tunneling matrix 
element, a\ a are the photon creation and annihilation 
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operators of the single- mode cavity with frequency uj, g 
is the qubit-cavity couphng constant, and ak{k — x, y, z) 
are the Pauh matrices. To facihtate the study, we write 
the Hamiltonian in the matrix form in units oifi = to = \ 



H 



a^a + g (at + a) — 



A 

2 



A 
2 



(2) 

To remove the hnear terms in at(a) a operators, we per- 
form the following two Bogoliubov transformations 



A = a + g,B = a — g. 



(3) 



In Bogoliubov operators A{B), the matrix element Hu 
{H22) can be reduced to the free particle number oper- 
ators A {B'^ B) plus a constant, which is very helpful 
for the further study. 

Different from the previous ansatz that the Hamilto- 
nian is expressed in the two operators A and B at the 
same time[l^, we here use the single operator one by 
one. First, in terms of operator A, the Hamiltonian can 
be written as 



„ , A'^A - a 



v4t A - 2g\A^ + A) + /3 ) ' 



where 



a = 5' + |, /3 = 3.92 + 1. 
The wavefunction is then proposed as 



(5) 



where e„ and /„ are the expansion coefficients, is 
called extended coherent state with the following prop- 
erties 

„„ _ (^')"|„)^.(^i±|):i„,^, 



10), 



e 2 



(7) 



The Schrodinger equation gives 



A 



n=0 



n=0 



n=0 



Left multiplying a ('tt-I gives 



A 



{m - a- E)em^ -jfm, 



(8) 



{m + P - E) fm-2g{m + 1) /„,+i - 2.g/„i_i = ye,„. 

(9) 

So the two coefficients e„ and with the same index n 
are related with 



A 



2{m-a- E)' 
and the coefficient /„ can be defined recursively, 

mf„r = ri(m - l)/,„-l - frn-2, 

1 / A2 

n{m) - (m + /3-i?)-- 

2g \ 4 [m — a — E) 



(10) 

(11) 
,(12) 



with /o = 1 and /i = f^(0) 

Similarly, by the second operator B, we can have the 
Hamiltonian as 



H = 
where 



^ f B>B + 2g{B^ +B)+/3' 

B^B -a' ' ' 



A 
2 



(13) 



a'=g'--, /3' = 3.g2 



2' ' " 2 
The wavefunction can be also suggested in terms of B as 



_r Er=o(-l)"V^/n 



(14) 



Proceed as before, the two coefficients and satisfy 
-2.9 (/^-i + {m + 1) /;,+i) + ("^ + /?' - i?) fL~^<n = 0, 



then we have 



m — a' — E 



f 



(15) 



00 

-2.g ^ (yfifnV^- \n - 1)^ + Vn+lV^.f„ \n + 1)^) 

00 

Ej2^-fn\n}A 



rj'(TO) 



r2'(m - l)/^_i - /;„_2: 
1 



2.9 



(to -I- ;3' - S) 



4 (m - a' - i?) 



(16) 
,(17) 



n=0 



with = 1 and /{ = n'{0). 

If both wavefunctions ([S]) and (jl4p are the true eigen- 
function for the same eigenvalue i?, they should be in 
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principle only different by a complex constant r if this 
eigenvalue is not degenerate 

oo oo 

Y^V^\e^\n)j, = rY,{-ir^/^\rjn)^, (18) 

oo oo 

Y,^\u\n)j, = vY^i-ir^KW)^. (19) 



n=0 



n=0 



Left multiplying the original vacuum state (0| to the both 
side of the above equations yields 

oo oo 

Y,V^\er.mn)^ = r^(-l)"V;^/;(0||n)^,(20) 



n=0 



n=0 

oo 



^V^/„(0||n)^ = r^(-l)"V^e;(0||n)^,(21) 
where 

^/^\{{)\n)^ = (-l)"V^(0|?i)5 = e-9'/2^". (22) 

To eliminate the ratio constant r, we have the following 
relation 



oo oo 



oo oo 



n— n— n— 

with the help of Eqs. ([TU| and ([T5|) . we arrive at 
A/2 



E 



;/nff"E 



n — a — E ^-^ n — a! — E 

n=0 n=0 



oo oo 



n=0 n=0 



(23) 



If set fn=K^ , f^=K^ and E = x — , we can recover 
Braak's exact solution [l3| 

G^{x) = (y)2i?+(a;)i?"(x) - R+{x)R-{x) = 0, (24) 
where 



R^ix) = E^"(^)3"' 

oo 

R^{x) = E 



n=0 



n=0 



a; — 71 =F 



If e = 0, the above equation can be reduced to the fol- 
lowing zero-bias case obviouslvfis) 

°° / A/2 \ 

G^(^) = E/"(^) lT^h/"-0. (25) 

n=0 ^ ^ 

Therefore Braak's G-functions are completely re-derived 
in a very intuitionistic and concise way. 

The G-functions can be written [19j in terms of so- 
called Heun functions [13] ■ However, the zeros of this 
Heun functions can not be given analytically, the numer- 
ical technique in the search for the zeros is still needed, 
so the cut-off for the summation should be unavoidable 
in the practical evaluation. 



B. Comparisons and discussions 

In the above derivation, the key point is proportional- 
ity of the two wavefunctions ([5] ) and (ITil) with the same 
eigenvalue. Both Hilbert spaces in the two Bogoliubov 
operators are complete, if truncation is not done, the 
proportionality is justified naturally for non-degenerate 
states. On the other hand, the degenerate eigenstates 
are excluded in principle here. It naturally follows that 
the Juddian solutions [21], l22| , which eigenvalue is doubly 
degenerate, are exceptional to Braak's solution. Interest- 
ingly, we do not need the extra condition of analyticity 
of the eigenfunction in Bargmann representation. In ad- 
dition, the validity of present approach is independent 
of the parity symmetry. The parity symmetry would be 
contained self-consistently in the final G-functions if the 
system really has, e.g. e — 0. 

Based on two Bogoliubov operator A and B, three 
present authors and one collaborator have used the fol- 
lowing wavefunction to the Hamiltonian(IT]) to analyze the 
spectrum in the qubit-oscillator systems [c.f. Eq. (6) in 
Ref. Kh]] 



J2n=0 Cn I") 



(26) 



where N is the truanted number. The numerical exact 
diagonalization (ED) in the space of the two Bogoliubov 
operators have given the spectrum exactly. The coeffi- 
cients c„ and dn can be obtained also. 

It is interesting to link coefficients in wavefunction (p6)) 
and those in wavefunctions (O and ([T^ as 



"n'.en, 



although the former ones are obtained from ED and the 
later ones by the zeros of the G functions. It can be also 
confirmed numerically. For practical interest, there are 
perhaps no essential differences between our work and 
Braak's solution fT3|, except that the avenues to obtain 
the basically same coefficients are different. The latter is 
described in a mathematical way and is of more concep- 
tual interest. 

In the mathematical sense, we can not rule out the 
possibility that the ED give good results for small N and 
get worse for higher N before the practical evaluation, al- 
though we know empirically that it is generally not that 
case for large N. For the low order perturbation the- 
ory, it happens that the third-order perturbation theory 
would give worse results than the second-order one in 
some parameter regime for instance, it may be not that 
case in very high order perturbation theory. In the calcu- 
lation, we really find that relative difference between the 
exact ones, which are easily obtained to any desired accu- 
racy, and those for the cut off N decreases monotonically 
with increasing N, and the convergence can be arrived 
at easily. One may know that the Heun series converges 
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A=1.0 



O Even parity with M=28 
O Odd parity with M=29 



LU 2 




FIG. 1: (Color online) Spectrum of the one-photon quantum 
Rabi model as zeros of functions F{a) in Eq. (|27|) at reso- 
nance. The numerical solutions are also collected with solid 
lines. 



before numerical calculations, although the cut-off can 
not be avoided in the real calculation. We indeed did not 
unfold the mathematics behind this formalism previously 
and focused on the analysis of the experimental spectrum 
data at that time[lo|. 

Braak's G-functions exhibit a very compact form in 
power series, which motivate us to reshape our previous 
work. By tunable extended bosonic coherent states, the 
QRM can be rnapped to a polynomial equation with a 
single variable |24j. We can also write this polynomial in 
power series as the following more concise form for large 
truanted number M 



F{a) 



M 

E 

J=0 



(27) 



where a is the key tunable variable we seek, and co- 
efficients are also related recursively with the following 
scheme 



, A, , , 



A 



2 ^ j\ 



(28) 



initiated from cq = 1.0 and ci = 0, because the coef- 
ficients with the two highest indices M and M — 1 are 
negligible small due to the required convergence and can 
thus be omitted. The zeros of the above function F{a) 
can also give the exact eigenvalue through 



+ A 



(29) 



where ± denotes the parity. The results are shown in 
Fig. m In Eq. (20) of the end of Ref. we have 



demonstrated that the wavefunction is equivalent to the 
expansion in the Fock space of displaced operators with 
tunable displacement. 

It is very interesting to note that zeros of the both func- 
tions defined through different power series in Eqs. ( [^5)) 
and (|27p can give the exact eigenvalues. In our practical 
evaluation, it is not more difficult to find the zeros for 
the function in Eq. (P7| than those in Eq. ( E5)) . because 
the poles at a; = n emerging in the latter are not present 
in the former. The key difference between Eqs. f B5t and 
(l27l) is that the former is well-defined without restriction 
and the latter is well-defined with build-in truncation. 



It is implied in the viewpoint that the QRM might 
have been solved exactly with an analytical closed-form 
solution in Ref. [l^. Nevertheless, whether Braak's ex- 
act solution could be called closed-form is subtle and 
therefore still controversial in our opinion. The so-called 
Heun functions can be basically called closed-form be- 
cause they are well defined, although much more com- 
plicated than e.g. the hypergeometric functions. But 
the eigenvalues are given by the zeros of the Heun func- 
tions, which can not be obtained without truncation in 
the power series. As shown in wavefunctions ([S] ) and 
([T4)). the expansion can not be closed naturally like in 
the JC model. It is generally accepted that the QRM has 
no trivial closed-form solution like that in the JC model 
due to the counter-rotating terms. The QRM can only 
have closed-form solutions with a vanishing qubit tunnel- 
ing A = OIU,!!^. Perhaps the question of " closed- form" 
solutions is academic and not of real value. Braak's so- 
lution is interesting for integrability of the QRM. 



C. Tutorial for Bogoliubov operators approach 

The present approach within the Bogoliubov operators 
can be generally described as follows, which is helpful 
for the further applications. The main task is to find 
the corresponding Bogoliubov operators. Then, one can 
expand the wavefunctions in terms of each Bogoliubov 
operator respectively. Eliminating the ratio constant of 
these wavefunctions will give the transcendental func- 
tions, which would be defined through power series in 
model parameter dependent quantities with coefficients 
related recursively. Finally, zeros of these transcendental 
functions would give the eigenvalues exactly, where nu- 
merical solutions to the one- variable (or finite variables in 
other multi-level systems for example) nonlinear equation 
must be required. Although the power series are defined 
through an infinite summation formally, in the practical 
calculation, it should be truncated to a finite summation. 
Fortunately, the obtained transcendental function G{x) 
can be written in terms of so-called Heun functions, by 
which we know the convergence before the numerical so- 
lutions. The unavoidable "cut-ofF' in the summation of 
the G-functions in the practical calculations means that 
some states in the Hilbert space is not considered, ac- 
cording to the wavefunctions ([SJ and (ITi|) . even they con- 
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tribute negligible small, they still belong to the Hilbert 
space of Bogoliubov operators. 

The applications to the two-photon QRM is performed 
in the next section as a preliminary example. 



III. EXTENSIONS TO THE TWO-PHOTON 
QRM 

The two-photon QRM is an phenomenological model 
describing a three-level system interacting with two pho- 
tons. When the intermediate transition frequencies are 
strongly detuned from cavity frequency, after adiabati- 
cally eliminating the intermediate levels, one arrives at 
the effective Hamiltonian. This two-photon QRM has 
also been studied for a long time with the RWA(26j and 
beyond the RWA[17|. [H, [Ig. It may describe the two- 
photon processes occurring in Rubidium atoms [13 
quantum dots 29] . 

The Hamiltonian of two-photon QRM takes the follow- 
ing matrix form 



H 



A 

2 



A 
2 



(at)%, 



First, we perform Bogoliubov transformation 

6 = ua + va) , 6^ = ua) + va, 



(30) 



(31) 



to generate a new bosonic operator. Compared to the 
Hamiltonian, if set 



/?+ 1 



(32) 



with P 



2 ' V 2 ' 
, we have a simple quadratic form of 



one diagonal Hamiltonian matrix element 



Hii = a* a + g 



Similarly, we can introduce another operator 

c — ua — va) , — ua) — va, 



(33) 



which yields a simple quadratic form of the other diago- 
nal Hamiltonian matrix element 



22 = a'a- g 



V 



In terms of the Bogoliubov operator b, the Hamiltonian 
can be written as 



H 



^ 2 
H22 



(34) 



with 



H'22 - {u^ 



+ Aguv) 
[uv + g {u^ + v^)] \{b^Y + b^ 



+ 2guv + V . 



The wavefucntions is suggested as 

( T,n=o ^/^e„ 1 71) J 



E„=0 ^-fn 



(35) 



where 



), = |0)6 = ^ rT^ \^)b : (36) 



6|0), = 0, 



(37) 



|0)j, is simply the single- mode squeezed vacuum state, 
is thus called as extended squeezed state. The 
Schrbdinger equation gives 

n=0 



n=0 



E Vri\en I n) 



b ' 



n=0 



n=0 

-[uv + g {u^ + v')] [{b^f + b'] J2 ^-fn \n\ 

n=0 n=0 



Left multiplying f, (to| gives 



- [uv + g {u^ + v^)] [/„_2 + (m + 2) (m + 1) J^+a] 
+ [(u^ +v'^)m + v'^ + 2guv {2m + 1) - £;] 
A 

-ye™ = 0. 

Then we have build one-by-one relation for coefficient e,^ 
and fm 



m — v^ 771 
— T-n — T ~^ -tv 



frr 



(38) 



which will considerably simplify the problem. The reclu- 
sive relation is then obtained as 



(m + 2) (to + 1) f.m+2 = -/m-2 + 



nim) 



where 



uv + g (u^ + V^) 

(39) 



Jl(m) = {u^ + v"^) m + + 2guv {2m + I) , 
A2 



i - e] 



(40) 
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The Hamiltonian can also be expressed in terms of the 
other Bogohubov operator c 



H = 



TTl A 

^11 



(41) 



\ 2 u^+v^ ) 

with 



+ Iguv + V 



The wavefunction then can be expanded in the Fock 
space of the c operator as the following form 



I Wn\e„ \n) 



(42) 



where I = n for n = 2k and I ~ n + 1 for n ~ 2k + 1. 
Therefore only two values of = ±1 are possible. 

Similarly, the Schrodinger equation gives the following 
relations 



- [5 {v^ + u^) + uv] [(m + 2) (to + 1) + /; 

+ [v^ + (u^ + v^)m + 2guv {2m + I) - E] 

d t 
f 



771 — 2 J 

A 



The coefficient Cm and fm are related by 

A 



and the reclusive relation is 



E 



(m + 2) (to + 1)/;„+2 
with 



/m-2 



^'(to) 



(43) 



-f 



UV -\- g (v? + v'^) 

(44) 



ri'(TO) = + {u^ + v'^) m + 2guv {2m + 1) 



-E - 



A2 



( 7n — v'^ 



E 



(45) 



Note that the two sets of coefhcients in the two wavefuc- 
ntions has the same form. 

Similarly, the two wavefunctions with the same eigen- 
value should be in principle proportional with each other 
for the non-degenerate state 



b - 



En=0 ' 



Left multiplying (0| to the both equations gives 

ri=0 ri=0 

Y,V^\U{Q\n), = rY,^'V^\e'^{(}\n)^. 



(46) 



We always have 



where 



=2k 



iWn\{0\n)^ = Vn!(0|n)j = L^'° 



{2k)\{uvf^ {-^y 



2k Z^^ji(^k~j)r 



T° 



{2k + l)\v{uvY 



J- 

. „ {2j + 1 



2'^j!(-^)^' 



2' ^ (2j + l)!(A:-j)! 



(47) 



(48) 



■(49) 



for even and odd particle numbers in the Bogoliubov op- 
erators b and c respectively. Then we have 



Now the summation is separated into two series with 
even and odd number n. To eliminate the constant r, we 
have 

/ . / 2 \ J" n / 2 \ " 

with the use of Eqs. §^ and (gS]). Set /„ = and 

—X ~ —v^ — E {u^ + u^) , we finally have 



, fn 



1 ± 



A 



2{n 



Ly = 0, 



(51) 



where the coefficient /„ is initiated from /o = 1 (/i — 1) 
for the case of the even (odd) n in the recurrence scheme 
Eq. ([5^ . and ± denotes the parity. Thus, G- functions 
for the two-photon QRM have been obtained. The zeros 
of the G-functions give the exact eigenvalues, as shown 
in Fig. [21 It should be straightforward to extend to the 
biased two-photon QRM, which is not shown here. 



Travenec 



18j has used Braak's approach to solve this 
model, but seems that the G-function like in Eq. ((?T|) 
is not given. Their coefficients are entangled in the two 
coupled equations, which may prevent such a simple de- 
scription for the G-functions. In the present Eqs. ( [38)) 
and ([^5)1 , two coefficients are related one-by-one with the 
same index n, which facilitates the derivations. This is 
also the advantage of Bogoliubov operators, which result 
in free particle number operators. 

The exceptional solutions to the two-photon QRM 
have been studied by Emary et al [2^. With the G- 
function in the two-photon QRM at hand, we can also 
discuss the Juddian solution similar to the one-photon 
model[l^. The G function is also not analytic in x but 

has simple poles at x = 0, 1, 2 For special values of 

model parameters g, there are eigenvalues which do not 
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0.1 0.2 0.3 0.4 0.5 

g 

FIG. 2: (Color online) Spectrum of the two-photon quantum 
Rabi model as zeros of G-functions in Eq. (|5ip . The numeri- 
cal solutions are also collected with solid lines. 

correspond to zeros of Eq. ([FT]) : these are the excep- 
tional sohitions. All exceptional eigenvalues are given by 
the positions of the poles 

E= (^n+i^ V-^ (52) 

which is exactly the same as that in Ref. [l^l- The nec- 
essary and sufficient condition for the occurrence of the 
eigenvalue is /n(x = n) = 0, which provides a condition 
on the model parameters g and A. They occur when 
the pole of Gf^g{x) a,t x — n is lifted because its nu- 
merator in Eq. (|5ip vanishes. From Eq. (|50|) . we know 
the proportionality is justified only for the even or odd 
photonic number respectively. The Juddian solutions are 
corresponding to those states which are degenerate, and 
therefore are excluded within this proportionality, so the 
lever crossing points with the same even, or odd photonic 
numbers are corresponding to Juddian solutions. 



IV. SUMMARY 

In this paper, by using the extended coherent states, 
Braak's exact solution in the QRM is recovered explicitly 
in an alternative more physical way. A preliminary ex- 
tension to the two-photon QRM is also performed with 
the use of extended squeeze states. The corresponding G- 
functions with the similar form of Braak's G-function are 
derived, which has not been obtained before. Both model 
can be treated in an unified way by the expansion in the 
Fock space of the Bogoliubov operators. Further exten- 
sions to other more complicated systems, such as multi- 
level, even multi-mode spin-boson model, are straightfor- 
ward, although perhaps a little bit tedious sometimes. 

For multi-level s pin- boson model, such as the finite- 
sized Dicke model[30|, the quantum chaos has been 
discussed sH . We have expanded the wavefunction in 

+ 1 Bogoliubov operators for Dicke model with fi- 
nite A^ two-level atoms, and get numerically exact so- 
lutions previously [ll]. According to the above discus- 
sions and the link with Braak's solutions, the exact solv- 
ability is ensured without doubt in this system. The 
quasi-integrability and the quantum chaos in this sys- 
tem should be very interesting. On the other hand, the 
multi-mode QRM has been also realized experimentally 
in circuit QED systems @. The extensions to these sys- 
tems are in progressing. 
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